A dissipative scheme to approach the boundary of two-qubit entangled mixed states 
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We discuss the generation of states close to the boundary-family of maximally entangled mixed 
states as defined by the use of concurrence and linear entropy. The coupling of two qubits to 
a dissipation-affected bosonic mode is able to produce a bipartite state having, for all practical 
purposes, the entanglement and mixedness properties of one of such boundary states. We thoroughly 
study the effects that thermal and squeezed character of the bosonic mode have in such a process 
and we discuss tolerance to qubit phase-damping mechanisms. The non-demanding nature of the 
scheme makes it realizable in a matter-light based physical set-up, which we address in some details. 
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The interplay between entanglement and mixedness 
has long been recognized as a crucial point to tackle 
towards a full understanding of the peculiar way cor- 
relations of a non-classical nature settle in multipartite 
configurations of quantum systems Despite the im- 
pressive efforts produced along these directions and a few 
important progresses being accomplished, a satisfactory 
comprehension of such an important topic is still elu- 
sive. Such incompleteness is well evident in the current 
unavailabity of a unique and unambiguous entanglement 
measure for general states involving more than two par- 
ties and the striking difficulties related to the ordering 
of mixed entangled states under different entanglement 
measures [2|. It is therefore of paramount importance 
to continue the investigation along these lines. An im- 
portant contribution to the problem represented by the 
trade off between entanglement and mixedness has been 
given by the classification of bipartite states exhibiting 
the maximum obtainable amount of entanglement for a 
given degree of mixedness 0, 0]. The explicit form of 
these genuinely interesting states, dubbed maximally en- 
tangled mixed states (MEMS's), strongly depends on the 
chosen measure for entanglement and mixedness. 

Devising means of realising these states is of great 
importance. In a way, as some form of noise will in- 
evitably be present in a physical set-up, the availability 
of pure states could well be out of question. The in- 
terest is thus in achieving the maximum possible entan- 
glement from the mixed-state resource one has to deal 
with. An efficient state-purification procedure can then 
be applied to MEMS's as described in Ref. [Bj]. Work 
at all levels has been performed on the generation of 
MEMS's. Linear-optics settings involving parametric 
down-conversion processes have been used in order to 
experimentally explore quite a substantial region of the 
physically allowed entanglement-mixedness space, up to 
the MEMS boundary @,Tzl- Theoretical proposals have 
been put forward for the navigation in the plane of entan- 
gled mixed states [H, Q, involving multi- level atom-like 
objects either interacting with structured environments 
or following properly arranged unitary evolutions. In- 
terestingly, Cho and McKenzie have proposed a scheme 
for the generation of bipartite Werner states in a two- 



impurity Kondo model via the well-known Ruderman- 
Kittel-Kasuya-Yosida interaction [10]. This is interest- 
ing as, under proper choices of entanglement and mixed- 
ness quant ificators, Werner states are MEMS's [3[. Here, 
we address a simple scheme that allows a system of two 
qubits to approach the boundary of physically allowed 
entangled mixed states. The protocol is based on bi- 
ased spin-boson interaction under the influences of dis- 
sipation and phase-damping. Besides its simplicity, the 
scheme addressed here shows that dissipation is able to 
coherently lead the qubit system to a partially entangled 
mixed state which, nevertheless, is interestingly close to 
the MEMS boundary. We show that the availability of a 
cold enough environment, together with an asymmetric 
preparation of the qubit system is all we need in order 
to approach such boundary curve. The scheme is highly 
realistic, as it estimates and includes the effects of the 
most relevant entanglement-spoiling mechanisms and, as 
we discuss, holds the promises for a prompt experimental 
realization in set-ups of cavity as well as circuit-quantum 
electrodynamics (QED). 



The remainder of this paper is organized as follows. 
Sec. U briefly introduces MEMS's and their properties, 
besides discussing the main technical tools used through- 
out the paper. In order to fix the ideas, Sec. Ill Al gives 
an account of the unitary version of the scheme discussed 
here while the main part of our analysis is presented in 
Sec. Ill Bi There, we address the reduced dissipative dy- 
namics undergone by the qubit system and identify a 
state of closest proximity to an element of the MEMS 
family. Some technical details, unnecessary to the com- 
prehension of our main results, are presented in Appendix 
A. In Sees. Ill Cl and III Dl we give account of how a struc- 
tured environment as well as the introduction of phase- 
damping mechanism would affect our findings. Sec. IIIII 
describes in some details an experimental set-up that 
has the necessary features for the implementation of the 
physical mechanisms assessed here. Finally, Sec. [IV] pro- 
vides a summary of our findings and conclusions. 
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I. INTRODUCTION TO MAXIMALLY 
ENTANGLED MIXED STATES 

Here, we briefly remind the basic properties of MEMS's 
and their parameterization for a specific choice of entan- 
glement and mixedness measure. As we mentioned, a 
long-standing argument regards state-ordering induced 
by entanglement measures [4|. It has been seen that 
for different entanglement measures there are different 
possible parameterisations of MEMS's. To date, bipar- 
tite MEMS's have been found to be described by a one- 
parameter family of density matrices [3[. For the pur- 
poses of this paper, we shall restrict ourselves to the 
entanglement measure given by concurrence [ill ], which 
is defined in relation to the entanglement of formation. 
The latter quantifies the number of Bell states required 
to prepare a given state. Entanglement of formation de- 
pends monotonically on concurrence which, for a two- 
qubit (pure or mixed) state p can be defined as [ll|, [l2[ 

C = max[0, - V% - - >Ad (1) 

where Ai > Xj (j = 2,3,4) are the eigenvalues of p(<j y <g> 
(jy)p*((jy (g) (jy) and a y is the y— Pauli spin operator. We 
shall be measuring mixedness of a state via linear entropy 
(see for example Bose and Vedral in [![) 

S=^[l-TV(p 2 )]. (2) 

The necessary parameterisation for MEMS with regards 
to these measures was provided in [3| , and identifies two 
subclasses of density matrices 
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where p x (p 2 ) holds for r e [2/3,1] (r e [0,2/3]). In 
a C versus S plane, these states lie on the so-called 
MEMS boundary curve shown in Fig. 1 (b), where the 
upper (lower) portion of the solid curve, i.e. the part 
corresponding to C > 2/3 (C < 2/3), is for pi (p2)- 



II. APPROACHING MEMS BOUNDARY 
A. Unitary case 

We start introducing the coupling model considered 
throughout our work, in the idealized situation of a per- 
fectly unitary evolution. Usually, problems involving 
the interaction of two-level systems with a boson can 
be modelled using an effective model where a dipole- 
like spin operator (proportional to the a x Pauli oper- 
ator) couples to the electric (or the magnetic) part of 
a field. This is the case for neutral atoms or quantum 



(a) (b) 




FIG. 1: (Color online), (a) Sketch of the discussed thought 
experiment. A cavity accommodates two two-level systems 
(labelled 1 and 2), asymmetrically coupled to the cavity field. 
The latter leaks out of the cavity due to a finite quality fac- 
tor, (b) Concurrence versus linear entropy MEMS boundary 
(full line). For comparison, the dashed lines shows the curve 
corresponding to the two-qubit Werner state pw = p\&) ($| + 
(l-p)t/4 with p e [0,1] and |$) = (l/>/2)(|00) + |H))i2. 

dots coupled to optical fields. However, this descrip- 
tion holds also for a system consisting of a Cooper-pair 
box [in a superconducting-quantum-interference device 
(SQUID) configuration and in the charge regime [14(] in- 
tegrated into a planar stripline resonator [15[, a setting 
generally referred to as circuit- QED. At the charge degen- 
eracy point, an effective dipole moment operator for the 
SQUID can be written, whose amplitude is proportional 
to the excess charge in the SQUID island [15]. Here, in 
order to fix the ideas and introduce the general formalism 
employed throughout our study, we use language and ter- 
minology typical of cavity-QED and we refer explicitly to 
a scheme of atomic qubits interacting with the field of an 
optical cavity. In Sec. [1111 we assess the details of possible 
experimental implementations. 

We consider two qubits with ground and excited states 
|0) and |1), respectively. They are allocated into a single- 
mode cavity and have the same transition frequency uo Ql 
resonant with the frequency uof ~ uo Q of the cavity field. 
This is described by the bosonic annihilation (creation) 
operator a {a)). Within the dipole-coupling interaction 
assumed here and using the rotating-wave approxima- 
tion [16] in interaction picture with respect to the free 
Hamiltonian of the system, the coupling reduces to (we 
assume units such that H = 1 throughout the paper) 

2 

W/ = 5Zft(^at + a+a), (4) 

where <t+ = (^jY = |l)j(0| is the qubit raising opera- 
tor and gj is the coupling strength of qubit j with the 
field. Fig. [1] (a) shows a sketch of the idealized situa- 
tion considered here. We assume that the qubit system 
is initally prepared in state |'0(O)) 12 = |01) 12 , while the 
cavity field is in |0) b . Eq. (jlj) commutes with the operator 
counting the total number of excitations within the sys- 
tem, so that the corresponding dynamics can be studied 
within finite-dimension subspaces. By considering the 
single-excitation subspace (consistently with our initial- 
state assumption) , it is straighforward to see that solving 
the time-dependent Schrodinger equation is equivalent to 
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finding the solution to the following set of coupled linear 
differential equations 



id t a(t) =gix(t), id t P(t) ■ 



(5) 



where we have used the decomposition a(t) |100) 12b + 
(3{t) |010) 126 + x(t) |001) 126 for the state of the whole sys- 
tem at time t (with |«| 2 + |/?| 2 + |x| 2 = !)■ By introducing 
the coupling ratio A = ^2/^1 and the dimensionless inter- 
action time r = git, the density matrix of qubits 1 and 
2 obtained by tracing out the field's state is found to be 



/X 2 (r) 



Qu(r) 



\ 










(6) 



with x(r ) = A sin(r yTTA 2 ") / vTTA 2 ", 0(t) = [1 + 
A 2 cos( rvTT A 2 ")]/(l + A 2 ) and a(r) = -A(l - 
cos(rvl + A^)(l + A 2 ). It is worth noting that in tracing 
out the field the observed dynamics are no longer evolv- 
ing unitarily with respect to r. Concurrence and linear 
entropy of g u (r) can be easily calculated. Their behavior 
is shown in Figs. [2] against r and A. By inspection, it is 
clear that concurrence and mixedness are maximized at 
small, non-unit values of the coupling ratio, when a sub- 
stantial amount of entanglement can be found in qubit 
states with S <C 1. Of course, mixedness of the state 
here arise in virtue of the loss of information over the 
field state. Shown in a C — S plane, these features make 
it evident that density matrices belonging to the bound- 
ary of physically meaningful states can be generated for a 
proper choice of r and A. Interestingly, despite the rather 
large temporal range we have considered (r G [0,40]), 
the elements of the MEMS family p2 cannot be produced 
by this scheme, which priviledges highly-entangled states 
(up to C = 1) of various mixedness. Differently from 
what was found in Ref. |8|, our states do not "track" the 
behavior of the MEMS boundary but follow a dynamical 
touch-and-go pattern with respect to p\. The similar- 
ity between g u (r) and p\ (for certain values of r, A and 
r) can be clearly understood by relying on the spectral 




FIG. 2: (Color online), (a) Concurrence of the produced 
q u {t) against the rescaled interaction time r = git and the 
coupling ratio A. (b) Mixedness of the state against git and 
A. 




FIG. 3: (Color online). Concurrence C versus linear entropy 
S in p u (r) and comparison with the MEMS boundary. Each 
continuous (blue) curve corresponds to an open trajectory of 
the two-qubit state having AG [0,2] (which increases at steps 
of 0.05) and the curvilinear abscissa is r G [0,40]. 



decomposition of these states: at the "touching" point 
between the open trajectories associated with two-qubits 
states shown in Fig. [3] and the MEMS curve, the eigenval- 
ues (eigenvectors) of g u are identical (locally equivalent) 
to those of pi [18]. This is never the case for states g u (r) 
and p2- 



B. Approaching MEMS boundary: Open-system 
dynamics 



Although promising because of the possibility of gen- 
erating a vast range of boundary states, the scheme 
described so far may be far from being realistic, in 
some set-ups. A unitary description is hardly retained 
for the whole range of time necessary in the MEMS- 
approaching mechanism when systems of quantum dots 
(neutral atoms) embedded in a cavity are considered, for 
instance. Any realistic physical set-up will imply the con- 
sideration of a finite rate of amplitude (phase) damp- 
ing affecting the system at hand. In particular, for the 
specific instance considered in this work, energy leak- 
age outside the cavity due to finite resonator quality fac- 
tor should be quantitatively included in our calculations. 
This is extremely important within the context of our in- 
vestigation, especially in virtue of the special role played 
by the excitation-conservation rules highlighted above. 
The introduction of an energy-dissipation mechanism, as 
it is the case with a leaky cavity, breaks such conserva- 
tion law forcing us to study the dynamics of the system 
in the whole Hilbert space, in principle. The effects of 
such differences should be carefully quantified. 

This is precisely what we do here, where we replace the 
Schrodinger equation at the basis of our study so far, with 
a dissipative master equation for the state of the qubits- 
field system g c (t) reading (unless otherwise specified, we 
use again the notation involving g\^ and t) 

d t g c (t) = -i[Hi,Q c (t)] +£[Q c (t)] = (4 +£)[Qc(t)] (7) 
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with the Liouville superoperator 

C[g c (t)} = 7 (n + l)(2ag c (t)a) - {a) a, g c (t)}) 
+ 7n(2a f ^ c (t)a - {aa 1 ", £ c (£)}). 



(8) 



Here, 27 is the energy dissipation rate from the cavity 
and n is the mean thermal occupation number of the bath 
with which the field is at equilibrium [19]. Throughout 
our treatment, we assume the bad cavity limit 7 ^> gi^- 
This implies that the cavity field response to the bath is 
much faster than that associated with its interaction with 
the qubits. In turn, this means that we can neglect the 
backaction of the photons emitted by the qubits therefore 
validating the Born approximation. Moreover, the bad- 
cavity limit also allows for Markovian treatment of the 
open dynamics as it is equivalent to a short qubit memory 
time. In analogy to what has been done in the unitary 
case, we want to trace out the field degrees of freedom, 
which can be systematically done by using projection- 
operator techniques [20[. However, in order to bypass 
the technicalities involved in such an approach, here we 
follow a simpler and yet rigorous way which leads to the 
same results as the projection-operator technique [2l|. 

In the bad cavity limit, it is reasonable to assume that 
the cavity field remains in the thermal steady state deter- 
mined by £, which we label p ss . We remind that within 
the bad cavity limit we are necessarily assuming weak 
coupling of the qubits with the field with respect to the 
field dissipation rate. In what follows, the validity of our 
results holds in this limit. Upon trace over the field state 
we get the following formal solution for the evolution of 
the qubits alone 

d t Q(t) = fTx^itue^-^tuW) ®p ss ]}dt f . (9) 
Jo 

Moreover, we can easily find that 



£>u(dQc) = 7(2n + l)dg c - 2jng c a, 
Cu(QcQ) = 27(n+ l)dg c - 7(2n + l)g c d 



(10) 



with g c = g (g) p ss . Upon iteration of these relations 
and explicit evaluation of the commutators involved in 
Eq. ([9]), a lengthy but straightforward calculation leads 
to the reduced master equation for systems 1 and 2 

d *e(t) = £ g M(n+ma-g(t)a+ - {a\al , g(t)}) 
j=i 7 

+ n(2a+g(t)aj-{atai,g(t)})} 
2 



fhfte 

■n(2a+g(t)a- -{&-&+, g(t)})}. 



£ {(n+l)(2a-g(t)a+ - {a+a~, g(t)}) 



(11) 



Eq. (pTj) is the starting point of our analysis. It can be 
solved by projecting it onto states of the two-qubit com- 
putational basis { 1 00) , |01) , 1 10) , 1 11) }i2 ? in a way so as 



to study the Bloch-like differential equations for the den- 
sity matrix elements g^i — 12(1 j\g\kl) 12 (i, j, fc, Z = 0, 1). 
The explicit form of such equations is given Appendix A. 
Here, we discuss the results achieved by solving them. 

Our initial state is again |01) 12 and the density matrix 
is found to have the general form 



/0oooo(t) 

00101 (*) Qouo(t) 

00110 (*) 01010 (*) 

V0 01111 (i) 



• (12) 



The presence of the non-zero piin(t) element is a first 
significant difference with respect to Eq. (j6]). In fact, 
while it is easy to see that in the unitary case concurrence 
was simply determined by the off-diagonal density matrix 
elements, here C depends critically on a delicate trade off 
between populations and coherences. Explicitly 



C 



-2g 0110 (t) - 2 v / £ooooO)£imO), (13) 



where we have used the fact that £0110 (i) < for any 
choice of the parameters involved and at any time. We 
now re-introduce the coupling ratio A and the dimension- 
less time r and analyze the behavior of C and S corre- 
sponding to state ([12]) . This is done in Figs. [4} where 
A G [0, 2] with r G [0, 100] are taken and concurrence 
(linear entropy) is studied for specific values of n. 

Immediately, one recognizes that C is maximized at 
A < 1, although the actual value is sensitive to n. Quite 
expectedly, in virtue of the features characterizing the 
unitary case, n — corresponds to the maximum (min- 
imum) of C (S) evaluated over the stationary state of 
the qubit system. This is reached already for r 30. A 
numerical inspection reveals that X opt ~ 0.8 corresponds 
to the largest possible value of concurrence of a state 
having S < 0.7. In the C — S plane, these results are 
summarized in Fig. [5] (a), where curves associated with 
increasing values of A and n are compared to the MEMS 
boundary. The (dashed) curve corresponding to X opt and 
n = is clearly highlighted. The tip of this curve corre- 
sponds to a state which is extremely close to the MEMS 




FIG. 4: (Color online). We show the behavior of linear en- 
tropy S (first row of plots) and concurrence C (second row) 
against A and r = git for n — 0, 0.4 and 0.8, in going from 
leftmost to rightmost plot in each row, respectively. 
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FIG. 5: (Color online). C — S plot for g(r) as determined 
by a dissipative dynamics with n > and j/gi — 10. The 
dashed curve embodies the open trajectory for the optimal 
state corresponding to n = and A = X op t- (b) Decay of 
entanglement and increase of mixedness in the optimal state 
for n that goes from to 1 (in steps of 0.1, from top to bottom 
open trajectory). 



boundary. However, a few remarks are in order: differ- 
ently from the unitary case, the dissipative scheme is able 
to approach family p 2l which belongs to the lower part 
of the MEMS boundary curve. Dissipation thus depletes 
the general properties of the generated states. Although 
very close to the boundary, q(t) never touches it. How- 
ever, the scheme does not require any time control. In 
fact, the state of closest distance from MEMS is achieved 
as a steady state: for r > 30 the properties of g do not 
change and the curves in Fig. [5] (a) do not fold back as 
time grows. 

The analytic form of the optimal state at any r can be 
found from the explicit solution of Eqs. (A-l) to (A-9). 

£oooo (t) = 7^ j £0101 (r) - 



G 2 



G 4 



£1010 (t) = — A , £0110 (r) = - V £ioio(t)£oioi(t) 



G 4 



(14) 



with £im (r) = and G = \J g\ + g 2 = #iVT + A 2 . It 
is interesting to notice that such an optimal state en- 
joys the same features as q u (t), i.e. the absence of 
populations of the |11) 12 state, which effectively leaves 
the two-qubit state in a one-excitation subspace. As 
for the unitary case, C is determined simply by the off- 
diagonal elements of the density matrix. We can quan- 
tify exactly how close we get to the p 2 class by using 
state fidelity F(pa,Pb) [12], which gives an estimate 
of the similarity between two density matrices pA and 
pB- For identical (orthogonal) states, F = 1 {F = 0). 
Here, we use the "amplitude" version of F defined by 
F(pa,Pb) = Tr^ ^J~pApB\fpA^ which is not too sensi- 
tive to very small changes in the optimal set of param- 
eters. By comparing p2 and the density matrix having 
elements given by Eqs. {HJ), we immediately recognize 
that F(g,p 2 ) cannot be close to 1 as the position of the 
coherence terms in the two density matrices do not cor- 
respond. Moreover, while r G [0,2/3] in p 2 guarantees 
positive coherences, we have that £ono(^) < Vr. And 



yet, the properties of the two states are rather close, as 
we have already commmented. However, it is straigthfor- 
ward to see that, by means of a bit and phase flip on qubit 
2, we get g p 2 . Needless to say, such local unitaries 
leave C or S invariant. The behavior of F(g, p 2 ) against 
r is shown in Fig. [6] (a) for increasing values of r in the 
range valid for p 2 . At r ~ 2/3, 7 = 10#i and A = 0.8 we 
find F(g,p 2 ) > 99.4%. This choice of parameters define 
the state (for r > 40) 



/0.398 



V 





0.362 
-0.295 






-0.295 
0.24 







0/ 



(15) 



which has concurrence and linear entropy values of C = 
0.589, S = 0.639. Fig. E (b) and Fig. 1 (b) show the 
effects induced of the thermal nature of the field on such 
an optimal state. In particular, state fidelity decreases 
quite rapidly towards ~ 0.7. Despite such a relatively 
large asymptotic value, this implies that the correspond- 
ing state is quite far from a MEMS, as seen in Fig. [5] (b) 
too. In fact, F = 0.7 is the state fidelity between the 
maximally mixed state 1/4 and p 2 . 



C. Irrelevance of the use of squeezed light for 
MEMS boundary approach 

Based on the results of Refs. [H, Q, where physical 
or effective squeezed fields where beneficial to the task 
of MEMS's generation, one may now wonder whether 
coupling the cavity field to a squeezed bath can improve 
the performances of the protocol addressed here. In this 
section we briefly show that, for the specific instance here 
studied, this is not the case at all. 

In order to provide a quantitative study to this sit- 
uation, we have derived the reduced master equation 
for two qubits immersed in the field of a cavity cou- 
pled to a broadband squeezed- vacuum characterized by 



(b) 




10 20 30 40 50 



-n,N 



FIG. 6: (Color online), (a) State fidelity F(g,p2) against the 
dimensionless interaction time r = git for A = X op t and r £ 
[0, 2/3] increasing at steps of 2/30 (going from bottom- to top- 
curve), (b) The dashed (solid) curve shows F(g,p2) against 
the parameter n (N) characterizing the thermal (squeezed) 
bath for r c± 2/3, git = 100 and A = X op t- The horizontal 
line shows the fidelity between p2 (with the above choice for 
r) and the maximally mixed state 1/4. 
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the parameters N and M < y/N(N + 1) related to the 
squeezing parameter s (the equality sign holding for ideal 
squeezing, in which case we have TV = sinh 2 ^). This 
can efficiently be done by writing Eq. (fTTj) for n = 
and replacing raising and lowering qubit opera tors w ith 
their "squeezed" versions = y/Na^ + y/ M/NcrJ 

and tj = y/M/Ncr+ + VN&J. The resulting master 
equation, depending explicitly on TV, M and the effective 
coupling rates #|/7 and #1(72/7, is identical to the one 
obtained adapting the approach described in [23[ which 
passes through the use of a squeezed and dissipative pic- 
ture. For M = 0, the thermal cavity- field case is found 
back, while the case of both M = N = corresponds 
to the optimal vacuum-field situation that has been ex- 
tensively discussed above. Bloch-like equations can be 
derived and an analysis analogous to the one described 
in Sec. Ill Bl performed. However, one recognizes that the 
introduction of a squeezed bath in the system considered 
in our study does not help the MEMS- approaching task. 
In fact, the situation described in Figs. [H seems to op- 
timize both the entanglement and mixedness within the 
two-qubit system. Without entering into the details of 
the analysis that has been performed, we simply mention 
that a numerical research of the point of closest approach 
to the MEMS boundary leads to the results shown by the 
solid curve in Fig.[6](b). These reveal that the largest fi- 
delity is achieved at N = (for an ideally squeezed field). 
The slower decay of fidelity in the squeezed case might 
be related to the protected steady state that is typically 
achieved when this sort of structured baths is considered. 



D. Inclusion of Phase Damping 

Until now we have only considered dissipation induced 
by losses in the cavity field due to its coupling with an ex- 
ternal bath. This resulted in effective damping terms in 
the reduced two-qubit master equation. However, spoil- 
ing mechanisms intrinsic to the qubit system, such as 
spontaneous emission and phase-damping, may originate 
important effects to be taken into account. Typically, 
these kick in with different time scales and their rel- 
ative weight is a system-dependent issue. However, a 
few general considerations can be made here. Sponta- 
neous emission from the excited state of each qubit oc- 
curring at rate j q <C gj <C 7 (in order for the bad cavity 
limit to hold) is formally accounted for by introducing in 
the two-qubit reduced master equation terms identical to 
the first line of Eq. (jTTJ) (for n = and replacing gj/j 
with j q ). This results in an effective modification of the 
spontaneous emission rate which becomes j q (l + Cj) with 
cj — g 2 jl (7^7) the (finite) cooperativity parameter corre- 
sponding to qubit j. Depending on the relative ratios of 
the parameters entering the problem, this simply biases 
the dynamics towards individual spontaneous emission 
processes of the two qubits, against the effective qubit- 
qubit interaction, therefore resulting in less entanglement 
and more mixedness of the stationary state. On the other 



hand, the inclusion of phase-damping processes in the dy- 
namics encompassed by Eq. ([TT]) needs to be commented. 

As the closest state to MEMS occurs for n = 0, we 
only need to consider this case. Phase damping is added 
by introducing 

C pd [g(t)] = -T[a z ,[a z ,g(t)]} (16) 

in Eq. ([7j). Notice that as C p d does not involve field op- 
erators, it remains unaffected by the procedure used to 
obtain the reduced qubit dynamics. Therefore, term ([T6]) 
is straightforwardly added to Eq. ([TT]) where we set n = 0. 
By means of suitably modified Bloch equations, we are 
able to quantify the effects that phase damping has on 
entanglement generation. Fig. [71 shows such effects on 
our optimal solution. We see that the inclusion of this 
term has a drastic effect on the state of closest proximity 
to MEMS. Even for the very small value T/g 1 w 0.001, 
we already see a marked drop in the amount of entan- 
glement our state possesses. In fact, Fig. [71 (b) shows 
that by T « 0.003 we have lost all the entanglement. 
Thus, a setting realizing the discussed scheme on effective 
MEMS's approach should keep phase-damping effects un- 
der control in order to mantain the protocol efficient. 

This result is interesting also under another viewpoint. 
In fact, it can be seen that despite the excellent proximity 
between p2 and our optimal state, the resilience proper- 
ties of the two states against phase damping processes 
are quite different: p2 as a family of states is incredibly 
robust to the action of phase-damping channels acting 
on the qubits. As shown by the dot-dashed part of the 
boundary curve in Fig. [71 (a), such a channel maps a state 
belonging to p2 into another one within the same fam- 
ily, so that the lower part of the MEMS curves basically 
folds on itself [24]. On the other hand, the optimal state 
g is affected in the rather different way addressed above, 
despite it is less than 1% away, in terms of state fidelity, 




FIG. 7: (Color online), (a) Effect of phase damping on the 
open trajectories corresponding to the states with optimal 
value X pt- The dashed curve shows the trajectory corre- 
sponding to T = 0, while the lower curves show an increasing 
dimensionless phase-damping rate T/gi which goes from 10 -3 
to 0.1 in steps of 10 -3 . The lower (dot-dashed) part of the 
MEMS boundary curve corresponds to MEMS p2 affected by 
independent phase damping channels, (b) Decrease in con- 
currence C against an increasing phase damping rate T/g± for 
the optimal state with A = X op t and git = 100. 
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from such a MEMS. Such a behavior asks for the design 
of suitable strategies for the protection of the proprties 
of g, at the optimal point, against processes bringing it 
away from the MEMS boundary. This is currently the 
topic of our ongoing investigations. 



III. PRACTICAL CONSIDERATION 

Here, we provide some details about an experimental 
setting we would like to suggest as a potential candidate 
for the successful implementation of the scheme we have 
discussed so far. It is somehow obvious that a cavity- 
QED set-up of two neutral atoms into an optical cavity 
as well as single-electron-charged quantum dots into a 
semiconductor microcavity would be well suited for the 
test of our proposal. However, here we decide to explic- 
itly approach a circuit- QED set-up, which we believe is 
the most promising scenario for our purposes in virtue of 
its properties of easy manipulability of stationary qubits 
and the physical features of the cavity and two-level sys- 
tems involved. 

In this scenario, each qubit is encoded into a standing- 
still superconducting qubit embodied by a SQUID work- 
ing in the charge regime at the degeneracy point (to 
wash away, to first order in the single- Cooper pair charge 
2e, the detrimental effect of low- frequency noise in- 
duced by background impurities) [3, El . Alternatively, 
one can use the recently proposed transmon qubits, a 
charge-phase qubit that results from a modification of 
a Cooper-pair box, which achieve umprecedented pro- 
tection against 1// noise. The transition energy of each 
superconducting qubit can be adjusted through an exter- 
nal, in situ magnetic flux that modulates the Josephson 
energy of the SQUID [14] in such a way that the qubit 
can be easily put in the strong resonant or dispersive 
regime with the field. This tuning ability is at the basis 
of the experimentally demonstrated non-demolition mea- 
surement of the qubit state through spectroscopic resolu- 
tion of the field's frequency-pulling effect [l5j |. The qubits 
are integrated, via conventional optical lithography, in a 
full- wave/half- wave on-chip coplanar waveguide split by 
input /output capacitances at tens of millimeters apart. 
The capacitors couple the cavity to input /output lines 
for the injection/leakage of the electromagnetic signal. 

The cavity resonance frequency is in the range of ~ 5 
GHz. With this frequency and an operating temperature 
< lOOmK, n is as small as 0.06, which allows for the 
vacuum- field treatment discussed in detail in this paper. 
Multiple-qubits have been experimentally allocated into 
a single coplanar waveguide, in a way so as to implement 
two-qubit information transfer via cavity-filed bus [26[ 
and, more recently, a circuit-QED version of the Tavis- 
Cummings model [25]. The biased coupling of qubits 1 
and 2 can be achieved by embedding them at slight asym- 
metric locations with respect to a voltage antinode of the 
sustained field mode. In principle, the stripline is a quasi- 
unidimensional structure with a very small transversal 



dimension that reduces the effective volume of the cavity 
field and enhances the coupling rate with the qubit. This, 
together with the effective dipole moment of the SQUID 
qubit (~ 2 x 10 4 eao) gives rise to gi/oof 0.2 — 2%. 
The energy damping time of the stripline can be as long 
as ~ l^s, which in principle allows for a long coherent 
dynamics within the cavity lifetime (experimental evi- 
dences put the qubit damping rate in the range of 2/is) 
and thus the implementation of the unitary version of 
the MEMS-approaching scheme. A detailed derivation 
of the qubit-stripline coupling Hamiltonian and the re- 
sulting coupling strength can be found in Paternostro et 
at. [15[. As the cavity into which the qubits operate is 
cut by interrupting the coplanar waveguide with the in- 
put/output capacitors, the resonator quality factor can 
be electrically tuned from 10 2 to 10 6 . This would in prin- 
ciple allow the realization of bad-cavity conditions where 
gi <C 7. Therefore, such a set-up is able to probe both 
the regimes studied here. Finally, we mention that for an 
individual-qubit it is T _1 ~ 2/is, which should put the 
phase-damping rate in a range of weak effect onto the op- 
timal state properties. Very recently, the ability to per- 
form complete state tomography of two transmon qubits 
has been experimentally demonstrated [2?}, opening up 
the possibility for preparation, evolution and characteri- 
zation of the target state. 

For the sake of completeness, we just mention that 
a different scheme for the preparation of MEMS using 
a pre-arranged off-line entangled resource has been sug- 
gested, both in cavity- and circuit-QED, in Ref. [9|. We 
refer to that work for an extensive account of the details 
necessary for such the step. 



IV. CONCLUDING REMARKS 

In this paper we have studied how dissipation can be 
used in a way so as to engineer a two-qubit state whose 
properties are fairly close to those of a MEMS. We found 
that concurrence as large as C « 0.59 can be set in a 
state having more than 99% fidelity with p2. This was 
achieved by introducing dissipation in a system than, in 
the noiseless unitary case, is able to span quite a large 
portion of the MEMS boundary. Although the noise- 
affected version of the protocol sees a strong reduction in 
its ability to produce genuine MEMS, yet a rather con- 
siderable possibility of navigation in the C — S plane is 
preserved, despite the explicit consideration of dissipa- 
tion. The relevance of our study is thus twofold. On one 
hand, it can be seen as the promising demonstration that 
a simple protocol for quantum state engineering of entan- 
gled mixed state can be designed for dissipation-affected 
settings. On the other hand, ours is the rigorous and 
complete quantitative assessment of the ability of a res- 
onant and bias spin-boson coupling to generate MEMS 
under the effects of relevant sources of noise. As such, 
we believe our investigation provides valuable and inter- 
esting information for the experimental groups interested 



in quantum state engineering in cavity and circuit- QED, 
where our study would find a natural and significant im- 
plementation. In perspective, it would be interesting to 
exploit the apparatus put forward here in order to ad- 
dress whether the system we have considered is able to 
create boundary entangled mixed states for more than two 
qubits, an issue which would require a considerable deal 
of theoretical work and represent a stimulating challenge. 
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APPENDIX 

As discussed in the main body of the paper, we can 
solve Eq. (11) by projecting onto the two-qubit compu- 
tation basis. In doing so we arrive at a set of differential 
equations that define the time-behavior of the density 
matrix. The amount of calculation required is greatly re- 
duced by exploiting the Hermitianity and normalization 
of the density matrix, thus reducing the set of relevant 
equations to 
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These equations, together with the normalization con- 
straint pun = 1 - (poooo + P0101 + P1010) and tne de- 
coupled equation c^poon = — ^~(2n + l)poon? allow for 
the solution of the dynamical problem discussed through- 
out the paper. 
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